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ABSTRACT 

A population of black holes (BHs) at high redshifts (z <; 6) that contributes significantly to the ionization of 
the intergalactic medium (IGM) would be accompanied by the copious production of hard ( <; 10 keV) X-ray 
photons. The resulting hard X-ray background would redshift and be observed as a present-day soft X-ray 
background (SXB). Under the hypothesis that BHs are the main producers of reionizing photons in the high- 
redshift universe, we calculate their contribution to the present-day SXB. Our results, when compared to the 
unresolved component of the SXB in the range 0.5-2 keV, suggest that accreting BHs (be it luminous quasars 
or their lower-mass "miniquasar" counterparts) did not dominate reionization. Distant miniquasars that produce 
enough X-rays to only partially ionize the IGM to a level of at most x e ~ 50% are still allowed, but could be 
severely constrained by improved determinations of the unresolved component of the SXB. 

Subject headings: cosmology: theory - quasars: general - intergalactic medium 
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1. INTRODUCTION 

The recent discovery of the Gunn-Peterson (GP) troughs in 
the spectra of z > 6 quasars in the Sloan Digital Sky Survey 
(SDSS; White et al. 2003; Wyithe & Loeb 2004) suggests that 
the end of the reionization process occurs at a redshift near 
z ~ 6. At this epoch, the ionizing sources drive a strong evolu- 
tion of the neutral fraction of the intergalactic medium (IGM) 
from values near unity down to x H i ~ 10" 3 (e.g., Fan et al. 

2002) . On the other hand, the high electron scattering optical 
depth, r e = 0. 17 ± 0.04, measured recently by the Wilkinson Mi- 
crowave Anisotropy Probe (WMAP) experiment (Spergel et al. 

2003) suggests that ionizing sources were abundant at a much 
higher redshift, z ^ 15. These data imply that the reionization 
process is extended and complex, and is probably driven by 
more than one population of ionizing sources (see, e.g., Haiman 
2003 for a post- WMAP review). 

The exact nature of these ionizing sources remains unknown. 
Natural candidates to account for the onset of reionization at 
z ~ 15 are massive, metal-free stars that form in the shallow 
potential wells of the first collapsed dark matter halos (Wyithe 
& Loeb 2003a; Cen 2003a; Haiman & Holder 2003). The com- 
pletion of reionization at z ~ 6 could then be accounted for by 
a normal population of less massive stars that form from the 
metal-enriched gas in larger dark matter halos present at z ~ 6. 

The most natural alternative cause for reionization is the ion- 
izing radiation produced by gas accretion onto an early popula- 
tion of black holes ("miniquasars"; see Haiman & Loeb 1998, 
Wyithe & Loeb 2003c, Bromm & Loeb 2003). The ioniz- 
ing emissivity of the known population of quasars diminishes 
rapidly beyond z <; 3, and bright quasars are unlikely to con- 
tribute significantly to the ionizing background at z <; 5 (Shapiro, 
Giroux & Babul 1994; Haiman, Abel & Madau 2001; Wyithe 
& Loeb 2003a). However, if low-luminosity, yet undetected 
miniquasars are present in large numbers, they could dominate 
the total ionizing background at z ~ 6 (Haiman & Loeb 1998). 
Recent work, motivated by the WMAP results, has emphasized 
the potential significant contribution to the ionizing background 
at the earliest epochs (z ~ 15) from accretion onto the seeds of 
would-be supermassive black holes (Madau et al. 2003; Ricotti 



& Ostriker 2003). The soft X-rays emitted by these sources can 
partially ionize the IGM early on (Oh 2001 ; Venkatesan & Shull 
2001). 

A population of miniquasars at z ^ 6 would be accompanied 
by the presence of an early X-ray background. Since the IGM is 
optically thick to photons with energies E below E m . dx = 1 .8[(1 + 
z)/15)] 5 xKi keV, the soft X-rays with E <, E m:ix would be 
consumed by neutral hydrogen atoms and contribute to reion- 
ization. However, the background of harder X-rays would red- 
shift without absorption and would be observed as a present- 
day soft X-ray background (SXB). In this paper, we examine 
the hypothesis that accreting BHs are the main producers of 
reionizing photons in the high-redshift universe, and calculate 
their contribution to the present-day SXB in this case. Our re- 
sults, when compared to the unresolved component of the SXB, 
suggest that accreting BHs cannot contribute significantly ei- 
ther to the completion of reionization at z > 6, or to the partial 
ionization of the IGM at z <; 15 to ionized fractions of x e J> 0.5. 

The outline of the paper is as follows: In § 2, we describe 
the method to calculate the SXB from quasars that contribute 
to reionizing the universe. In § 3, we critically discuss cur- 
rent X-ray observations, focusing on the unresolved fraction of 
the SXB that could be attributed to distant quasars. In § 4, we 
calculate the expected contribution to the SXB from hypothet- 
ical quasars and their lower-mass miniquasar counterparts that 
fully ionized the universe. In § 5, we repeat our analysis for 
a putative miniquasar population that partially ionizes the IGM 
at high redshifts. In § 6, we discuss how various simplifica- 
tions made in our analysis influence our final results. Finally, in 
§7 we summarize our results and the implications of this work. 
Throughout this paper, we adopt the background cosmological 
parameters as measured by the WMAP experiment, fl,„ = 0.27, 
n A = 0.73, n b = 0.044, and h = 0.71 (Spergel et al. 2003) and 
set the mass fraction of helium to Yn e = 0.24 ( e.g. Buries, Nol- 
lett, & Turner, 2001). In the rest of the paper 'accreting BHs' 
will refer to both quasars and their lower-mass "miniquasar" 
counterparts. 
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2. MODELING THE CONTRIBUTION TO THE SXB 

To calculate the contribution to the present-day SXB from 
accreting BHs that reionized the universe, we assume for sim- 
plicity that the accreting BHs form in a sudden burst at redshift 
z = Zq. The spectrum of the ionizing background is a crucial 
ingredient of the modeling, and depends on the type of accret- 
ing BH that is considered. The details of the assumed emission 
spectrum for each type are discussed in §4 and 5. In § 6, we ex- 
amine the dependence of our results on the emission spectrum. 

Assuming that the cumulative hard X-ray flux of the accret- 
ing BHs results in a present-day X-ray background F E , we can 
compute the number of H-ionizing photons present at redshift 
Zq per unit comoving volume: 

n =— dE— (1) 

c Jn± eV E 

l+Zg 

In the above equation, E and F E correspond to photon energies 
and intensities at z = 0, and F E has the units of ergs cm" 2 s" 1 
deg" 2 erg _1 . The photon energy E m . dx = 1.8[(1 +z)/15)] ' 5 keV 
reflects the value above which the IGM becomes optically thin 
across a Hubble length for a fully neutral IGM 1 . For most spec- 
tral slopes in the range we study here, we have verified that our 
results are insensitive to the choice of E mux . The exception is in 
§5, where we consider unusually hard spectra; in this case the 
exact value of E max does influence the results. A brief discus- 
sion of this will be given in §5. 

The total number of H-ionizations per unit comoving volume 
is approximately equal to the comoving density of H-ionizing 
photons n 7 with energy below £ max . This simple formulation 
ignores secondary ionizations by the photo-electrons produced 
by the more energetic photons (E > 100 eV). For spectra of the 
form F E oc E~ a with a > 0.7 at E > 13.6 eV, these secondary 
ionizations are subdominant compared to the ionizations by the 
UV-photons (Abel et al. 1997). More importantly, the num- 
ber of secondary ionizations is negligible for any spectrum if 
the ionized medium has a neutral fraction, xm <C 1 (since the 
electron energy is instead used to heat the gas; e.g. Shull & van 
Steenberg, 1985). For harder spectra and a partially ionized 
medium, secondary ionizations become important. We will re- 
turn to this issue in § 5 below, where we discuss a miniquasar 
population that is assumed to only partially ionize the IGM. 

Since helium may be reionized later by the already known 
population of optical quasars at z ~ 3 (Sokasian, Abel & Hern- 
quist 2002; Wyithe & Loeb 2003a), we consider only hydrogen 
here. The presence of helium will not significantly change our 
results, as will be discussed in §6.2 below. 

The normalization of the present-day background F E is ob- 
tained by requiring that the quasar population would lead to 
reionization, 

« 7 = V n H.O (2) 

where n H ,o = Q h h 2 p all (\ - Yn e )/m p = 2.05 x 10" 7 cm" 3 is the 
number density of H-atoms at z = and rj is the number of ion- 
izations per hydrogen atom that are required to achieve reion- 
ization. 

The value of 77 is larger than unity because of the need to 
balance recombinations. The ratio between the Hubble time, 

1 More precisely, the upper limit should be taken as the energy below which 
photons would typically be absorbed and contribute to H-reionization at any 
redshift prior to z ~ 6 for a partially ionized IGM - rather than E max as defined 
above (see §5, eq. 10). In practice, the difference between the two definitions 
is small. 



which is taken to be fhub(z) = 2//" 1 (z)/3, and the recombina- 
tion time in a region at overdensity 5 that is partially ionized 
to xh+ = nH*/riH, where tih is the number density of hydrogen 
atoms (neutral and ionized) at any redshift z, is given by: 

Where x e = n e /riH (= xh+, because helium is not included), 
C = (n 2 )/(n} 2 which depends on the level of clumping in the 
high-redshift the IGM, and is expected to be > 10 (Haiman, 
Abel, & Madau 2001; Wyithe & Loeb 2003a; Ricotti & Os- 
triker 2003), and T is the gas temperature. This implies that 
recombinations cannot be ignored at high redshift. In fact, the 
number of recombinations per hydrogen atom between z = Zq 
and z = 7, N KC , can be estimated from 

f' (z=7) a„{T) n e n H+ 

Nrec = / dt (4) 

Jt(z=zg) H H+ 

Interestingly, for any combination of zq and x e yielding the 
same r e , this can be written as a function of r e as long as x e is 
constant with redshift for z > 7 and x e <~ 1 for < z < 7: 

where r e is the electron scattering optical depth between red- 
shifts z = and z = Zq. Note that x e <~ 1 for < z < 7 and x e = 
for z > 7 yields r e = 0.05. The recombination coefficient cch(T) 
is taken from Abel et al. (1997). To maintain a fixed ionization 
fraction we need r\ = N rec + 1 <~ 25 ionizations per hydrogen 
atom. In §6.5 we will show that the actual value of 77 depends 
on the exact topology of reionization, and therefore quite un- 
certain. The SXB intensity F E predicted below simply scales 
linearly with the exact value of 77. Given the uncertainty in 77 
we will adopt the fiducial value of 77 = 10, and scale our quoted 
fluxes in units of (77/IO) (see Table 2). 

It is important to emphasize that we are normalizing F E based 
on the total ionizing background, required to produce reion- 
ization. Our prediction for F E therefore does not require the 
knowledge of the luminosity function of the quasars. 



3. THE UNRESOLVED SXB AND ITS UNCERTAINTIES 

Since the unresolved fraction of the SXB is the key observ- 
able used to obtain the constraints in this paper, we first discuss 
its measurement and uncertainty. The soft X-ray band in the 
energy range 0.5-2.0 keV has been studied by Moretti et al. 
(2003, hereafter M03). We focus on this energy range since 
it provides the strongest constraints. The 0.5-2 keV band 
corresponds to hard X-rays at z > 6, which are well within 
the regime of photon energies to which the IGM is optically 
thin. M03 determines the intensity of the total SXB within 
this band to be / F E dE ~ EF E (1 keV) = 7.53 ± 0.35 x 10" 12 
ergs cm" 2 s" 1 deg" 2 , when combining 10 different measure- 
ments reported in the literature. M03 further include deep pen- 
cil beam surveys together with wide field shallow surveys and 
ultimately are able to identify discrete sources in the flux range 
2.44 x 10" 17 - 1.00 x 10" 11 ergs cm" 2 s" 1 within the soft band. 
They find that 94^% of the SXB is made up of discrete X- 
ray sources. The majority of this resolved component involves 
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point sources (88 ±7% of the total SXB), with extended sources 
amounting to 6% of the background. Barger et al. (2002, 2003) 
show that the bulk of these sources are at low redshifts, z < 4. 
In § 6.4 below we find that any z > 6 quasars that "hide" among 
the already resolved point sources can only make up to 1 % of 
the total SXB. 

The fraction of the SXB contributed by known discrete (point 
or extended) sources depends on the sensitivity limits of cur- 
rent observations. M03 find that when they extrapolate the an- 
alytic form of the log N-log S distribution (cumulative frac- 
tion of the total SXB in sources brighter than flux S) to S m j n ~ 
3 x 10" 18 ergs s" 1 cm" 2 , which is a factor of ~ 10 lower than 
current sensitivity limits, the entire SXB can be accounted for 
by discrete sources. The additional, still undetected faint X- 
rays sources could be either point-like or extended. 

Wu & Xue (2001) calculate the theoretically expected amount 
of X-ray emission in the soft band from thermal emission by 
gas in all clusters and groups in the universe using the ob- 
served L x -T relation and X-ray luminosity function (XLF) 
and find this should be 1.18 x 10" 12 ergs cm" 2 s" 1 deg" 2 . This 
is 16% of the total SXB quoted by M03. Approximately ~ 
40% of this contribution originates in galaxy groups; ignor- 
ing groups (where the XLF extrapolation becomes uncertain) 
provides a conservative estimate of the total expected contri- 
bution from thermal emission from extended sources, <~ 0.6 x 
1.18 x 10" 12 = 0.68 x 10" 12 ergs cm" 2 s" 1 deg" 2 , which is 9% 
of the total SXB. This suggests that the actual contribution to 
the SXB from extended sources is higher than claimed by M03, 
and that deeper X-ray observations will reveal new extended 
X-ray sources at flux levels below current sensitivity limits. To 
remain conservative we take the M03's estimate of 6% for the 
contribution of extended sources to the SXB, but keep in mind 
that the actual number is probably higher, leaving a smaller 
residual that can be attributed to distant accreting BHs. 

Recent numerical simulations by Dave et al. (2001) suggest 
that a fraction as large as 0.3-0.4 of all baryons might be in 
the warm/hot intergalactic medium (WHIM) with temperatures 
between 10 5 - 10 7 K, residing in diffuse large scale structures 
such as the filaments connecting virialized objects at z, = 0. The 
WHIM should emit some thermal radiation in the soft X-ray 
band, but has so far not been detected. In this paper, we ignore 
the possible contribution to the SXB from the WHIM. Note that 
any WHIM contribution would strengthen our results, since it 
would again leave a smaller residual of the SXB that can be 
attributed to distant accreting BHs. 

Soltan (2003) pointed out that point sources can never make 
up the full SXB, because the X-ray photons will be Thompson 
scattered by free electrons, yielding a diffuse component. This 
component is estimated to be <~ 1 .0- 1 .7% of the total SXB. Al- 
though this is a small fraction, it makes a non-negligible contri- 
bution to the unresolved SXB, and we include it in our analysis 
below. 

A summary of the possible SXB contributions is given in 
Table 1. The Table quotes an intensity for a "mean" and a 
"maximum" unaccounted SXB flux. The mean unaccounted 
flux is the total flux minus the contributions from point sources, 
galaxy clusters, and Thompson scattered X-rays. The maxi- 
mum unaccounted flux is calculated by summing the lcr errors 
(3 rd column in Table 1) in the most conservative way: the total 
SXB is set to its maximum (+lcr) value, while all other com- 
ponents (point sources, extended sources and scattered X-rays) 
subtracted from this total are assumed to be at their minimum 



(—lcr) levels. We find the mean and maximum unaccounted 
fluxes to be 0.35 x 10" 12 ergs cm" 2 s" 1 deg" 2 and 1.23 x 10" 12 
ergs cm" 2 s" 1 deg" 2 , respectively. 

It should be noted that a similar analysis can be done for the 
hard X-Ray background (2.0-8.0 keV), but since the constraints 
in this band are considerably weaker, we do not discuss them 
here. 



Table 1 

Contributions to the Soft X-ray Background 



Source 



Flux 

(10" 12 ergs crrf 2 s" 



1 deg" 2 ) 



Total SXB 1 

Point Sources 1 

Diffuse Component 1 

Expected scattered background 2 

Unaccounted flux (mean) 

Unaccounted flux (max) 



7.53 
6.63 
0.45 
0.10 
0.35 
1.23 



0.35 
0.35 
0.25 
0.03 
0.60 



'From Moretti et al. 2003. 

2 Diffuse Thompson-scattered radiation from point sources (Soltan 2003) 



4. FULL REIONIZATION BY THE UV EMISSION OF 
ACCRETING BLACK HOLES 

4. 1 . Reionization by Quasars 

We follow the procedure outlined in §2 and calculate the 
present day SXB, Fe, for a burst of quasars at redshifts zq = 
6,10 and 20, assumed to fully reionize the universe. We adopt 
the shape of the fiducial quasar spectrum from Sazonov, Os- 
triker, & Sunyaev (2004, see their eqs. 8 and 14), which is 
then assumed to redshift to z = without absorption for hard X- 
rays. This broad-band spectral energy distribution of the aver- 
age quasar in the universe is computed based on a detailed theo- 
retical model that fits a range of data, including published com- 
posite spectra of quasars in the optical, UV and X-ray bands. 
The template spectrum scales approximately as Fe oc E~ 11 for 
E > 13.6 eV, but becomes much shallower beyond E = 2 keV, 
where it scales as Fe oc E~° 25 . The main reason for this 'kink' 
in the spectrum around E = 2 keV is to accommodate the rise 
in the observed cosmic X-ray background in the energy range 
1-30 keV. 

The resulting flux density is plotted in the range 0.1-10 keV 
in Figure 1, assuming i] = 10, with the upper curve correspond- 
ing to zq = 6, and the lower curve to zq = 20. The EFg curves 
are redshifted and suppressed by a factor of (1 +zq) relative to 
the rest-frame template spectrum (which peaks at E = 60 keV). 
As a result, we obtain the largest SXB for the quasar burst at 
the low end of the redshift range, zq = 6. 

The two straight lines plotted in the range 0.5-2.0 keV denote 
the mean and maximum unaccounted flux of Table 1 . The un- 
accounted flux is assigned a power-law spectrum with a slope 
matching the spectral slope of quasars in the range 0.5-2.0 
keV. The figure shows that at all three redshifts, the SXB in- 
tensity significantly exceeds the level of the mean unaccounted 
flux (by a factor of ~ 4-6), and is above our inferred maximum 
level. To compare directly with the values quoted in Table 1 we 
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calculate Fqso = Ji^ST F E dE. For z = 6, 10, and 20, we find 
F QS0 = 2.17, 1.84, and 1.37 x 10 _12 (7//10) ergs cm" 2 s" 1 deg" 2 , 
respectively. For any redshift, the maximum unaccounted flux 
is saturated. 




E (keV) 



FIG. 1 . — Predicted soft X-ray background due to a burst of quasars produc- 
ing r) = 10 ionizing photons per hydrogen atom and reionizing the universe at 
z = 6, 10 and 20 (first three curves from the top). The template spectrum for 
these quasars is taken from Sazonov, Ostriker, & Sunyaev (2004). The two 
straight lines plotted in the range 0.5-2.0 keV are the maximum and mean un- 
accounted flux in the observed soft X-ray background (Table 1). A power-law 
spectrum is assigned to this flux, with a slope matching the spectral slope of 
quasars. 



4.2. Reionization by Miniquasars 

In the previous section, we addressed the scenario in which 
the IGM is fully ionized by luminous quasars. An alternative 
scenario is that in which the ionizing background is dominated 
by accretion onto the much less massive seed BHs that later 
grow to be the massive BHs powering luminous quasars. The 
main difference, for our purposes, is the typical spectrum of 
the ionizing sources in the two scenarios. Our fiducial tem- 
plate from Sazonov, Ostriker, & Sunyaev (2004) that we used 
above is appropriate for luminous quasars at lower redshifts 
(0 < z < 5), powered by supermassive black holes with masses 
of Mbh ~ 10 6-9 M Q . The spectra of miniquasars, with BHs 
whose masses are in the range M bh ~ 10 2 " 4 M (Madau et al, 
2003) are likely to be harder. 

To obtain the spectrum of a typical miniquasar, we follow 
the approach of Madau et al. (2003) and adopt a template that 
consists of two components. One component originates from a 
multi color accretion disk (MCD), which describes the accre- 
tion disk as series of blackbody annuli with different temper- 
atures, and results in a spectrum that scales as Fe ,mcd oc E 1 ' 3 
up to kT ma!l = IkeV (M^/Mq)' 1 ^. The second component is 
a simple power-law Fe.pl oc E~ l for E > 13.6 eV. The precise 
origin of this power-law emission is not fully understood, but 
is likely to be a combination of Bremsstrahlung, synchrotron, 
and inverse Compton emission by a non-thermal population of 
electrons. 



This power-law + MCD composite model adequately describes 
some ultraluminous X-ray sources (ULXs), that are claimed to 
be associated with an intermediate mass black hole of ~ 10 3 
M Q (Miller et al. 2003). In these spectra the flux from the 
MCD, Fmcd , is comparable to the flux in the power-law com- 
ponent, Fpl. Miller et al. (2003) find the ratio of Fp^ to the total 
flux to be 0.33 and 0.63 for the two ULX sources NGC 13 13 X- 
1 and X-2 respectively. Motivated by these observations, here 

we take Fmcd = J ™ 6keV ^k,MCD^, ^pl = J 2 keV F E ,phdE, 
and define the ratio = Fmcd/^pl- In this paper, we adopt 
\£ = 1, and vary it between 0.5 and 2.0 as an additional source 
of uncertainty 2 . 

An additional source of uncertainty is the mass of the typical 
black hole, Mbh- Madau et al. (2003) show the mass function 
of accreting intermediate mass black holes at four different red- 
shifts. The median mass at z = 22 is Mbh ~ 100 M Q , increasing 
to ~ 1000 Mq at z = 13. To approximately mimic this increase, 
we adopt log(M b h) = (4.1- 0. lz) ± 0.5. The error in the black 
hole mass is incorporated into the error in the final contribution 
to the SXB. 

Applying the analysis outlined in § 2 with this modified spec- 
tral shape yields F m Q S o = J^skeV ' F E dE = 3 - 3 ± L1 > 2.3 ±0.8 
and 1.4 ±0.4 x 10 _I2 (?7/10) ergs cm" 2 s" 1 deg" 2 for z = 6, 10 
and 20, respectively. The dominant contribution to our quoted 
error is the uncertainty in "f, the relative importance of the flux 
in the power law component to the flux from the MCD compo- 
nent. Despite the differences in the assumed spectra for quasars 
and miniquasars, their contribution to the current SXB, under 
the assumption that they fully ionized the IGM, is comparable. 
For any redshift, we find that miniquasars saturate the mean un- 
accounted flux of the SXB. Only for z = 20 the maximum unac- 
counted flux of the SXB is not saturated within the uncertainty 
of our model. 

5. PARTIAL "PREIONIZATION" BY X-RAYS FROM 
MINIQUASARS 

We have so far omitted secondary ionizations from our anal- 
ysis. Secondary ionizations are collisional ionizations by the 
photo-electron that is knocked off the hydrogen atoms when 
ionized by energetic X-ray photons. Shull & van Steenberg 
(1985) quantify the fraction <p(xn+,E) of the energy of these 
electrons that is used for H-ionizations. They find that for 
energies ^100 eV, 4>( x h+,E) is independent of energy, and 
4>(xh+,E) ~ 1/3 for a neutral medium (xh+ = 0) although it 
drops rapidly for Xh* > 0.1. This justifies for our statement 
in § 2 that secondary ionizations do not contribute significantly 
to full reionization (as we further show in more detail at the end 
of this section). However, secondary ionizations can become 
significant in a partially ionized medium, especially for harder 
spectra. 

The large electron scattering optical depth WMAP result can 
be explained with a partially ionized IGM, followed by full 
reionization at z < 6.3 (e.g. Wyithe & Loeb 2003b, 2004; Cen 
2003b; Haiman & Holder 2003). For example, reionization 
up to an ionized fraction x e = 0.3,0.5. and 0.7 at zq = 34,24, 
and 20, respectively, with x e = const down to z ~ 7 and x e ~ 1 
at < z < 7 yields r e = 0.17 (Ricotti & Ostriker 2003). In 

2 The definition of Fpl is the same as in Miller et al. (2003), whereas our 
definition of Fmcd is not exactly the same because our expression for Fe mcd is 
only an approximation to the actual spectrum (see Fig 3. in Shakura & Sunyaev 
1973). The error this introduces is well within the uncertainties we assign to 
the value of . 
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these preionization scenarios, secondary ionizations are impor- 
tant because the quasar spectrum is significantly harder and 
consists primarily of X-rays. Since helium is ignored, x e = x H +. 
In §6.2 we demonstrate that including helium does not change 
the final results by more than ~ 15%. 

X-rays may have some advantages in partially reionizing the 
universe (Oh 2001; Venkatesan, Giroux, & Shull 2001; Ri- 
cotti & Ostriker 2003). First, for low ionization states of a gas 
( x h+ ^ 0.1), a single X-ray photon can produce multiple sec- 
ondary ionizations. Also, the escape fractions of X-rays from 
the local high density environments of sources into the lower 
density regions, is high: / esc ~ 1 . In these low density regions, 
the clumping factor C is expected to be <~ 1 . 

This implies that, contrary to the case of full ionization by 
UV-photons, which requires 77 > 10 photons per baryon (§2), 
partial ionization to a fraction x H + by X-ray photons requires 
only r] = x H +(l +N rec )/(1 + J\f ssx ), where 7V se c is defined to be 
the mean number of secondary ionizations caused by a single 
photon (given below in eq. [9]). In the preionization models of 
Ricotti & Ostriker (2003) a considerably harder template spec- 
trum is used, with F E oc E 4 at E < 100 eV, and oc E~ l at E > 100 
eV, as appropriate for highly obscured sources. A physical mo- 
tivation for the break at this energy is that the column density 
of isothermal gas in hydrostatic equilibrium in a typical virial- 
ized halo at z = 10-20 is JVhi ~ 10 20 " 21 cm -2 (Dijkstra et al. 
2004). Photons with E > 150 eV (N ffl /10 20 cm" 2 ) 1 / 3 can es- 
cape unimpeded from a halo, whereas the halo is optically thick 
to less energetic photons. 

In order to repeat the analysis of §2 with the modified spec- 
trum and the secondary ionizations included, we first compute 
the mean number of secondary ionizations per photon, Af sec - 
The total number of ionizations per unit comoving volume to 
go from electron fraction x e to x e + dx e is: 



dn Hion = (1 +N lec ) dx e n H ,o 



(6) 



The number of hydrogen ionizations per comoving volume 
is given by: 



"Hion 



= n 7 ( 



1 + 



1 1+z 4i 
n-v 13.6 eV c 



Jm e v 



dE F E x 



x0[£(l+z)-13.6eV,x e ]) =n 7 [l + (Mec)(x e )] (7) 

In this equation, n 7 is the comoving number density of H- 
ionizing photons (eq. 1). As mentioned above, the more en- 
ergetic photons of energy E will produce a secondary photo- 
electron of energy £-13.6 eV, which will spend a fraction 
4>(E- 13.6 eV,x e ) of its energy on H-ionizations. The fitting 
formula we used for <P(E,xh+) in our numerical integration is 
described in the Appendix. Note that we can only write x e here 
because it is equal to x H + (see §6.2 for the effect of helium). The 
quantity (Af sec ) (x e ) denotes the mean number of secondary ion- 
izations per photon over the entire energy range at electron frac- 
tion x e . The above equation simply states that each H-ionizing 
photon produces one primary and (A/" se c) (x e ) secondary ioniza- 
tions, when the electron fraction in the gas is x e . Note that the 
primary ionizations are actually driven by helium, but as argued 
in §6.2 this only slightly changes our final results. 

Combining equations (6) and (7), we require for the number 
of ionizing photons per comoving volume 



(1+jVrec) n H fi dx e 
1 + (Kcc)(Xe) 



The total number of ionizing photons per comoving volume, 
« 7 , tot, to change the ionized fraction from to x e can be ob- 
tained by integrating equation (8) from to x e . 



I7, tot — 



(1+ Nrec) x e n Hfi 
1+Mec 



(9) 



in which the number Af sec is the mean number of secondary ion- 
izations per photon over the whole energy and electron fraction 
range. This is the number we referred to earlier in this section. 

To show the relative importance of secondary ionizations, we 
list the values of Af sec in Table 2. Note that these numbers are 
slightly higher than those in Madau et al. (2003). The main 
reason is that our number is mean number of secondary ioniza- 
tions per photon over the whole electron fraction range up to x e , 
whereas their number is given at the specific value x e . 

Given A/" S ec, we set rj = x e (l+N lec )/(l+Af sec ) an d calculate the 
present-day SXB intensity for the three preionization scenarios 
with different x e and reionization redshift, as considered by Ri- 
cotti & Ostriker (2003). The fluxes we find are listed in Table 2. 
As the table shows, models within which X-rays are assumed 
to partially ionize the IGM to x e < 0.5 do not saturate the mean 
unaccounted SXB. The preionization scenario with x e = 0.7 has 
a similar contribution to the SXB as the models that aim for full 
reionization and saturates the maximum unaccounted flux. 

As mentioned in §2 the results are dependent on the exact 
value of £ max , being the energy below which photons would 
typically be absorbed and contribute to H-reionization at any 
redshift prior to z ~ 6, for a partially ionized IGM. The optical 
depth of a partially ionized IGM between z = Zq and z = 6 for a 
photon of energy E emitted at z = Zq is, 

T(E): 



f a dz 

■■ / -rn[n HI a H (E')+n HeI a He (E')], 



(8) 



(10) 

where E' = E(l +z)/(l +zq) and n Hl (n He i) is the proper number 
density of neutral hydrogen (helium) atoms. Although helium 
has been ignored in this section, it is included in this formula 
since it dominates the contribution to t(E) for X-ray photons. 
We discuss the implications of helium in more detail in §6.2. 
The cross-sections are taken from Verner et al. (1996). 

The value of E m - dx was obtained by setting, T(E msa ) = 1 . There 
is some arbitrariness in this definition, and to check how it af- 
fects our results we repeated our analysis using E miix:d \ t defined 
by r(£ , max a i t ) = 0.1. This roughly corresponds to doubling £ max 
and brings down our calculated contributions to the SXB by 
< 15%. Since this choice of £ max ,ait is very unrealistic (photons 
of energy £ max ,ait are only absorbed ~ 10% of the time) we be- 
lieve that the actual choice of Zs max is influencing our results by 
less than <~ 15%. 

Using the model for secondary ionizations, we can calculate 
their contribution for the full reionization cases described in §4. 
We find that for the quasar case (§4.1), A/" sec £ 0.04, which is a 
negligible contribution. For the miniquasars, with harder spec- 
tra, A/^ec S 0.10, which would come down to a ~ 10% reduc- 
tion in the calculated contribution to the SXB. Since the errors 
in these contributions were of the order of ~ 30%, ignoring the 
contributions from secondary ionizations is justified. 

6. DISCUSSION 

6.1 . Constraints on Spectrum 

As emphasized above, the spectral shape of the putative typ- 
ical high-z accreting BH is uncertain; the existing templates, 
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motivated by lower-redshift sources, can be considered merely 
as guides. In this section, we adopt a different strategy and ig- 
nore any physically motivated features in the spectrum. Instead, 
we adopt a simple power-law shape, F E oc E~ a for E > 1 3 .6 e V. 
We then consider a range of values for a, and also vary rj (eq. 2) 
to calculate the contribution to the SXB. Figure 2 shows which 
combinations of a and rj value (Eq. 2) are allowed by the mean 
and maximum unaccounted flux in the SXB (see Table 1). The 
solid (dotted) lines denote those models with combinations of 
a and rj such that the predicted contribution to the SXB is equal 
to the mean (maximum) unaccounted flux. The redshift of the 
burst of accreting BH formation in this case is assumed to be 
ZQ = 10. 

This figure shows that for -q = 10, which was used in §4.1 
and §4.2, any power-law shallower than aw 1.4 (1.2) will sat- 
urate the mean (maximum) unaccounted flux. This figure also 
demonstrates that when a = 1, the mean (or maximum) unac- 
counted flux is saturated when rj is greater than 1 (or 4). 




FIG. 2. — Constraints on the number of ionizing photons per H atom, r;, 
(eq. 2) and the power-law index for the slope of the ionizing background, a, 
based on the intensity of the present-day SXB. The quasars are assumed to 
form at zn = 10 an d have a power-law spectrum, Fe oc E~ a for E > 13.6 eV. 
The curves bracket the allowed parameter space for the mean (solid line) or 
maximum (dotted) unaccounted flux in the SXB. 



6.2. The Influence of Helium 

So far our analysis has ignored the presence of helium for 
simplicity. In this section we will demonstrate how the presence 
of helium reduces the contribution of accreting black holes to 
the SXB, albeit at a level < 15%. 

The photoionization cross-section for neutral (and singly ion- 
ized) helium, an e (and cxHe+) is significantly larger in the X- 
rays than that for hydrogen, a#. For example, for E = 100 eV 
CHe/cH ~ 20 and for greater energies rapidly increases to ~ 30 
after which this value stays constant. Since the number density 
of helium atoms is <~ 8% that of hydrogen, photoionization of 
helium will dominate those of hydrogen by a factor of <~ 2-3. 
The primary photo-electrons created will cause secondary ion- 
izations of mainly hydrogen according to the description in §5. 



For simplicity, let us first neglect recombinations and assume 
that only helium is photoionized once, after which the gener- 
ated photo-electrons spend 36% of their energy on hydrogen 
ionizations and 5% on helium ionizations (Shull & van Steen- 
berg 1985). With our assumed spectrum and taking zq = 34 
we obtain 8.3 secondary ionizations of hydrogen and 1 . 1 of he- 
lium per photoionization of helium. This implies that for every 
ionized helium atom, 8.3/2.1 = 4 hydrogen atoms are ionized. 
Since the hydrogen atoms outnumber the helium atoms by a 
factor of <~ 12, x He + rises faster than xh+ by a factor of <~ 3 for 
a neutral medium. 

This ratio is reduced when direct photoionization of hydro- 
gen is taken into account. For photons more energetic than 
100 eV, roughly 1 hydrogen atom is photoionized per 2 pho- 
toionized helium atoms. This results in xn e + rising faster than 
x H + by a factor of ~ 2.5. The number density of electrons is 
n e = n H + + n He + + 2n He ++. Initially n e = n H + + n He + = n H (x H + + 
0.08x// e +) ~ 1 .2hh* (since xh c +/xh+ <~ 2.5 and xn e ++ = 0). 

In fact, the abundance of He" 1 " 1 " is negligible even when he- 
lium starts to become significantly ionized, because the recom- 
bination time of He ++ is shorter than that for hydrogen by a 
factor of ~ 5. Venkatesan, Giroux, & Shull (2001) included 
helium in their models and found that xn e ++ ~ 0.01xf/ e +. When 
He ++ recombines, a UV-photon will be emitted that can ion- 
ize hydrogen and helium with roughly equal probability, since 
0He ~ 14<7h for a 55 eV photon. Similarly He" 1 " will recombine 
N lec times emitting a UV-photon, which is more likely to ion- 
ize hydrogen than helium by about a factor of 2. This implies 
that a fraction ~ 2N lec /3(l+N Tec ) of the photons that ionized 
helium will end up ionizing hydrogen. This will bring down x e 
to <~ \.\5xh*. Once He" 1 " 1 " recombinations become more abun- 
dant, the ratio x e /xn* will be reduced further, since the pro- 
duced UV-photon can ionize hydrogen. Therefore, the main 
effect of recombinations is to make x He + closer to x H + and to 
bring x e closer to ~ 1 . lx H + which were the underlying assump- 
tions in the models of Shull & van Steenberg (1985). This jus- 
tifies our use of their fitting formula <P(E,xh+). 

To estimate how the above effects influence our results in §5 
we replace x e in equation (7) by 0.87x e (because x e ~ 1 . 15xh+ as 
argued above) and also multiply the number of secondary ion- 
izations by a factor of 1 . 14 to account for secondary ionizations 
of helium occur at a level of 14% relative to that of hydrogen 
(Shull & van Steenberg 1985). This generally leads in a re- 
duction of at most ~ 15% to the contribution to the SXB. The 
'helium corrected' contributions are given in Table 2. Since 
the helium correction is small, the main conclusions will not 
depend on whether helium is included or not. 

The < 15% helium correction should also not be taken too 
seriously since the method outlined in §5 is highly simplified. 
We did not follow the temperature evolution of the gas, which 
is a strong function of the ionized fraction of the gas. For low 
ionized fraction x H + < 0. 1 , T is expected to be much lower than 
10 4 K, whereas it may become much higher once the gas is ion- 
ized; the value of T affects the recombination rates. Moreover, 
the incident spectrum at a certain location will be a function 
of time and distance to the nearest source, since an interven- 
ing column of neutral hydrogen and helium (which change as a 
function of time) will harden the spectrum. A proper treatment 
of this problem requires 3D numerical simulations with radia- 
tive transfer and goes beyond the scope of the present paper. 
Keeping in mind the simplifications inherent in our discussion, 
the helium correction is most likely within the uncertainty of 
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the model. 

In summary, our estimates of the present-day contribution 
to the SXB presented in §5 were slightly too high, because we 
ignored helium. We demonstrated that the amount we overesti- 
mated this contribution is less than a factor of ~ 1 . 15, which is 
most likely within model uncertainties. 

6.3. Effects of Continuous (Mini)Quasar Formation 

Our analysis has so far required that all quasars at z = Zq pro- 
duce r) ionizing photons per hydrogen atom, as the condition 
for achieving complete reionization. A more realistic model 
would include a continuous formation of accreting BHs over an 
extended redshift interval. There is no reason to assume that 
their formation is stopped once the number of photons required 
to ionize the IGM is produced. Our previous conclusions are 
therefore conservative, in the sense that the accreting BHs that 
form between z m in <^ z < Zq would represent an additional con- 
tribution to the unresolved SXB. Here z m ; n <~ 4 is the highest 
redshift at which faint quasars would be resolved in existing X- 
ray surveys (Barger et al. 2002, 2003); quasars at z < z m in are 
already included in the resolved fraction of the SXB. 

The presence of accreting BHs forming beyond those needed 
for reionization will increase the total quasar contribution to 
the SXB and tighten our constraints. Naively, the correction for 
this effect would increase with zq, since the range z m ; n <; z < 
Zq available for the continued formation of accreting BHs is 
longer. As was shown in Figure (1) and Table 2, the contribu- 
tion to the SXB from an ionizing background is higher at lower 
redshift at a fixed value of -q. Modeling the magnitude of the 
correction due to the extra accreting BHs is beyond the scope 
of this paper; we merely emphasize here that the presence of 
these extra sources would tighten our constraints. 

6.4. Known X-ray Point Sources 

As discussed in §3, a fraction of 88 ± 7% of the SXB is ac- 
counted for by known point sources. We next discuss whether 
some of these point sources may in fact be unidentified high 
redshift quasars. The expected X-ray flux from an accreting 
BH at redshift zq, shining at the Eddington Luminosity, for a 
black hole of mass Mbh is, 

-20 / M bh 



S QSO = 3.8x 10" 



V10 4 M, 



.( fx 
A 0.03 



10 



(l+zfl)-l 



9/4 



ergs sec 1 cm 2 



(11) 



where fx is the fraction of the total emitted energy in the (0.5 - 
2.0)(1 +Zq) keV band. The value f x = 0.03 corresponds to the 
template spectrum as given by Sazonov, Ostriker, & Sunyaev 
(2004). The last term is a fit to d 2 L {\0) / dl(z Q ), where cI l {zq) is 
the luminosity distance to z = Zq (Pen, 1999). The fit is accurate 
to <~ 4% over the range of 6 < z < 34. 

For reasonable choices of parameters, we do not expect any 
of the known point sources to be the low-mass miniquasars 
we considered above. However, massive black holes M bh > 
10 8 Mq with harder spectra could have detectable fluxes, and 
could be included in the sample of known faint X-ray point 
sources. The work by Barger et al. (2002, 2003) suggests that 
most of the known point sources in the Chandra Deep Field 
South (CDF-S) are at redshifts z < 4. Nevertheless, Koekemoer 
et al. (2004) find 7 faint sources in the CDF-S, with counter- 
parts in near IR JHK-bands with no flux in any of the bands 
up to 8500A. Koekemoer et al. (2004) speculate that these 7 



sources with extreme X-ray/optical flux ratios (EXOs), could 
be high-redshift quasars (6 < z < 10). However, the total back- 
ground flux from these seven EXOs is 8.4 x 10 -15 ergs s _1 cm -2 
in the 0.5-8.0 keV band, which corresponds to an SXB flux 
of 7.6 x 10 -14 ergs s _1 cm -2 deg -2 if such sources populate the 
whole sky. Therefore, we conclude that even if all of these 
sources turn out to be quasars at z > 6, their contribution to the 
SXB is still < 1% of the total, increasing the allowed SXB from 
the total high-redshift accreting BH population by < 25%. 

6.5. Reionization Topology and r\ 

We mentioned in §2 that the exact value of r\ depends on 
the topology of reionization. We saw that r\ is very different 
between a case in which reionization is mediated by UV and 
X-ray photons. In the first case, initially 77 > 10 because the 
UV-photons cannot escape the high density regions in which 
the effective clumping factor C (eq. 4) is boosted by small scale 
density fluctuations to values C > 10 (Haiman, Abel, Madau 
2001). After the filaments are cleared of their neutral hydrogen 
gas the UV-photons can escape freely into the low density re- 
gions, in which C <~ 1 . In this type of "inside-out" scenario, the 
clumping factor is a decreasing function of cosmic time, and r\ 
is large ( <; 10). In the second case, X-rays are able to escape 
the filaments unimpeded and immediately start partially ioniz- 
ing the low-density voids. Reionization may then be completed 
by stars producing UV photons. This scenario does not, in fact, 
require hard spectra, as long as the sources send most of their 
ionizing radiation directly into voids, rather than ionizing local 
dense regions. The clumping factor in this "outside-in" case is 
an increasing function of cosmic time, and 77 can be low (e.g 
Miralda-Escude et al. 2000), as long as most of the mass at 
z <~ 7 is still not highly ionized. Both scenarios are theoreti- 
cally plausible, and depend on the nature and the location of 
the first ionizing sources. Although it is known that in the local 
(z ~ 0) universe, galaxies lie along dense filaments, the exact 
distribution of the high redshift ionizing sources is uncertain 
and will influence -q. 

7. CONCLUSIONS 

Our main results, obtained in § 4.1 and § 4.2, are summa- 
rized in Table 2. The third column in this table denotes the per- 
centage of the total SXB contributed by accreting BHs (which 
we used to refer to both quasars and their lower-mass "mini- 
quasar" counterparts) for models in which high redshift accret- 
ing BHs contribute significantly to full or partial reionization. 
These percentages can be compared directly with the percent- 
ages given by M03. The second column gives the total inten- 
sity of the quasar populations, to be compared to the allowed 
range of the mean or maximum unresolved SXB flux of 0.35 x 
10 -12 ergs cm -2 s _1 deg -2 to 1.23 x 10 -12 ergs cm -2 s _1 deg -2 , 
respectively. 

As Table 2 shows, models in which z > 6 accreting BHs con- 
tribute to reionization overproduce the SXB. Pre-ionization by 
miniquasars requires fewer ionizing photons, both because only 
a fraction of the H atoms need to be ionized, and because hard 
X-rays can produce multiple secondary ionizations. We find 
that models in which X-rays are assumed to partially ionize the 
IGM up to x e <~ 0.5 at 6 < z < 20 are still allowed, but could 
be severely constrained by improved determinations of the un- 
resolved component of the SXB. 

We emphasize that our constraints derive from the total num- 
ber of ionizing photons that the population as a whole needs to 
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produce. Therefore, our conclusions depend mostly on the as- 
sumed spectral shape, and are independent of the details of the 
population, such as the luminosity function and its evolution 
with redshift. Future improvements in resolving the SXB, im- 
proving the limits on the unresolved component by a factor of 
a few, would place stringent constraints on the contribution of 
z ~ 15 accreting BHs to the scattering optical depth measured 
by WMAP . 
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APPENDIX 

We used a fitting formula for the fraction <fr(x, E) of the en- 
ergy of a fast photo-electron used for H-ionizations, which is a 
modified version of the formula provided by Shull & van Steen- 
berg (1985). Their work focusses mostly on energies E 3> 0.1 
keV, in which <fi(x,E) is independent of energy. The energy de- 
pendence at lower energies, 0.05 < E < 0.5 keV, can be seen 
most clearly in their Figure 3, to which the formula below pro- 
vides an accurate fit: 



a(x,E) = — arctan 

7T 



0(x, J B) = O.39[l-x - 4092 ^] 1 ' 7592 
/ E \/0.03 \ - 25 
V012keV J V^5~ J . 



(12) 



At high energies (E > 0.5 keV) this approaches the fitting 
formula provided by Shull & van Steenberg (1985). Note that 
the E = 28 eV curve is fitted less accurately for x > 0.2. The 
photons in this energy range however, are not affecting our anal- 
ysis in §5. 
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Table 2 

Unaccounted Soft X-ray Background (0.5-2.0 keV) vs. Quasar/Mini-Quasar Contribution 



Source 


Flux 

(10~ 12 ergs cm -2 s -1 deg~ 2 ) 


percentage (%) 
of total SXB 


A/sec 


Unaccounted flux (mean) 0.35 ± 0.60 


4.5 ± 8.0 




Unaccounted flux (max) 1.23 


16.3 






Quasars §4.1 








Units of Flux here (j?/10) 10~ 12 ergs cnT 2 s 


- 1 deg- 2 




zq = 6 


2.17 


29 




ZQ = 10 


1.84 


24 




zq = 20 


1.37 


18 




Mini-Quasars §4.2 


ZQ = 6 


3.3 ± 1.1 


44 ±14 




ZQ = 10 


2.3 ±0.7 


30±11 




zq=20 


1.4 ±0.4 


19±5 






Pre-ionization Scenarios §5 








Units of Flux here (^ffi^) 10~ 12 ergs cm" 


2 s- 1 deg- 2 




x e =0.3, zq 


= 34 0.22 


3 


2.44 


x e =0.5, zq 


= 24 0.72 


10 


1.50 


x e =0.7, zq 


= 20 1.5 


20 


0.97 


Helium correction applied to models above (see §6.2). 


x e =03, zq 


= 34 0.19 


3 


3.03 


x e =0.5, zq 


= 24 0.61 


8 


1.94 


X e =0.7, ZQ 


= 20 1.3 


17 


1.30 



